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Abstract 

By using the Lindblad theory for open quantum systems, an analytical ex- 
pression of the tunneling probability through an inverted parabola is obtained. 
This penetration probability depends on the environment coefficients. It is shown 
that the tunneling probability increases with the dissipation and the temperature 
CO ■ of the thermal bath. 
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Quantum tunneling with dissipation has been intensively investigated in the last two 
decades [1—12]. Very interesting is the discussion whether the dissipation suppresses or 
enhances the quantum tunneling. Caldeira and Leggett [H [2] concluded that dissipa- 
tion tends to suppress quantum tunneling. Using a different method, Schmid obtained 
similar results in Ref. |13j . Widom and Clark [H] considered a parabolic potential 
barrier and found that dissipation enhances tunneling. Bruinsma and Bak [15] also 
considered tunneling through a barrier and found that at zero temperature the tun- 
neling rate can be either increased or decreased by dissipation. Leggett [IE] considered 
tunneling in the presence of an arbitrary dissipation mechanism and found that, nor- 
mally, dissipation impedes tunneling, but he also found an anomalous case in which 
dissipation assists the tunneling process. Razavy [17] considered tunneling in a sym- 
metric double-well potential and concluded that dissipation can inhibit or suppress 
tunneling. Fujikawa et al. [18] also considered tunneling in a double-well potential 
and found an enhancement of tunneling. Harris [19] calculated the tunneling out of 
a metastable state in the presence of an environment at zero temperature and found 
that quantum tunneling is enhanced by dissipation. In [20], Yu considered the tun- 
neling problem in an Ohmic dissipative system with inverted harmonic potential and 



he showed that while the dissipation tends to suppress the tunneling, the Brownian 
motion tends to enhance it. In a series of papers [2TJ [22], Ankerhold, Grabert and 
Ingold have studied real time dynamics of a quantum system with a potential bar- 
rier coupled to a heat bath environment, employing the path integral approach. The 
conclusion drawn from these papers is that different authors have studied different 
problems using different methods. They obtained results which in many cases present 
an evident contradiction. 

In the present paper we study the tunneling in the presence of a dissipative en- 
vironment in the framework of the Lindblad theory for open quantum systems, based 
on completely positive dynamical semigroups [23-25]. We extend the work done in 
some previous papers [26j [27] . In [27] a similar problem was treated by using the path 
integral method and numerical calculations. Our study can be applied to problems of 
nuclear fragmentation, fission and fusion, considered as a tunneling process through 
the nuclear barrier defined in the space of collective coordinates, like charge and mass 
asymmetry or the distance between the fission fragments. 

For the inverted harmonic potential, the tunneling problem in the framework of 
the Lindblad theory can be solved exactly. In Sec. 2 we write the basic equations 
of the Lindblad theory for open quantum systems and give results for the coordinate 
and momentum expectation values and variances for the damped inverted harmonic 
oscillator. Then in Sec. 3 we consider the penetration of a Gaussian wave packet 
through the potential barrier and define the penetration probability. In Sec. 4 we 
analyze its dependence on various dimensionless parameters which enter the theory 
and show that the probability increases with the dissipation and the temperature of 
the thermal bath. A summary is given in Sec. 5. 

2 Quantum Markovian master equation in 
Lindblad theory 

The simplest dynamics for an open system which describes an irreversible process is 
a semigroup of transformations introducing a preferred direction in time [23-25]. In 
Lindblad's axiomatic formalism of introducing dissipation in quantum mechanics, the 
usual von Neumann-Liouville equation ruling the time evolution of closed quantum 
systems is replaced by the following quantum master equation for the density operator 
p(t) in the Schrodinger picture [21], which is the most general Markovian evolution 
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equation preserving the positivity, hermiticity and unit trace of p: 

^ = ~[H, P(t)} + i £(2^ - V?ty>(t) - p(t)V?yj). (1) 

Here if is the Hamiltonian operator of the system; Vj, Vj are bounded operators on 
the Hilbert space of the Hamiltonian and model the effect of the environment. We 
make the basic assumption that the general form ([I]) of the master equation with a 
bounded generator is also valid for an unbounded generator. 

As usual, we define the two possible environment operators V\ and V2, which are 
assumed linear in momentum p and coordinate q, as follows [28, 29J: 

Vj = a 3 p + b 3 q, j = 1, 2, (2) 

with a,j, bj complex numbers. The Hamiltonian H is chosen of the general form 

H = H + ^(qp + pq), H = ^-p 2 + U (q) , (3) 
where U(q) is the potential. With these choices and with the notations 

h h h 

fl w=n E \ a j\ 2 , D PP = 9 E N ; A><? = D ip = -9 Re £ a*bj,\ = -Im Y, a ?i' ( 4 ) 
z j=l,2 z 3=1,2 z 3=1,2 j=l,2 

where a* and b* denote the complex conjugate of ctj and bj, respectively, the master 
equation (jTJ takes the following form [28, 29J: 

dp % % % 

-j t = —AH ,p] - ^(A + p)[q,pp + pp] + ^(A-p)[p,pg + gp] 

-^[Q, [q,p]] ~ ^[P, \p,p}} + ff%, Ml + \P, hp}})- (5) 

Here the quantum diffusion coefficients D pp , D qq , D pq and the dissipation constant A 
satisfy the following fundamental constraints [2B1 29J: D pp > 0, Dqq > and 

DppDgq - D 2 pq > (6) 

In the particular case when the asymptotic state is a Gibbs state 

Pg(oo) = e kT /Tre *=t ; (7) 

the coefficients for a harmonic oscillator potential have the following form [28. [29]: 

A + p Tioj A — jj, h Tioj 

Drjr, = nrnio coth — — -, D aa = coth — — , D va = 0, (8) 

pp 2 2fcT' m 2 moo 2kT' pq K J 

3 



where T is the temperature of the thermal bath. The fundamental constraint ([6]) is 
satisfied only if A > [i. 

In the following we denote by oaa the dispersion (variance) of the operator A, i.e. 
oaa =< A 2 > — < A > 2 , where < A >= a a = Tr(pA) is the expectation value of the 
operator A and Tip = 1. By a AB = 1/2 < AB + BA > - < A >< B > we denote the 
correlation of the operators A and B. 

From the master equation (jSj) we obtain the following equations of motion for the 
expectation values and variances of coordinate and momentum: 

^- = -(X-p)a q (t) + ^a p (t), (9) 

^-<^>-(a + ^«) do) 

and, respectively 



da qq (t) 



2(A - n)a qq {t) + -a pq {t) + 2D qq , (11) 



dt K r ' WK ' m 

*± = _ 2(A + , )<rw(t) _ < ^ P + P f^M > +2 < f^M > , (() + 2Bw , (12) 

%W = - < d JM q > + < > „,(() + - 2Xa„(t) + 2D„. (13) 

at aq dq m 

For the harmonic oscillator with the potential U(q) = muj 2 q 2 /2, the solutions of 
these equations of motion are obtained in Refs. [28l [29] . In this paper we consider the 
tunneling through a potential barrier given by an inverted harmonic potential (inverted 
parabola) with 

U{q) = ~^fq 2 . (14) 

The Hamiltonian Hq in equation ([3]) with the potential (fl4l) can be regarded as the 
Hamiltonian of a harmonic oscillator with an imaginary frequency iu and the equa- 
tions of motion fl9l- fTl~3l) for this potential are formally obtained by performing the 
replacement u — > iu in the corresponding equations for the harmonic oscillator. These 
equations of motions can be solved by using the same method as in references [281 129] 
for the harmonic oscillator. Contrary to the situation of the harmonic oscillator, where 
we have two cases, overdamped and underdamped, for the inverted parabola such a 
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distinction does not exist. The solutions for the expectation values and variances of 
coordinate and momentum coincide formally with the solutions corresponding to the 
overdamped case of the harmonic oscillator. For the expectation value of the coordi- 
nate and momentum we obtain with v = y/u 2 + fi 2 : 

a (t) = e~ xt ((coshut+ - sinh vt)aJ0) + — sinh vtaJO)), (15) 

v mv 

2 

a p (t) = e~ A *( sinh uta q (0) + (cosh vt - - sinh ut)a p (0)). (16) 

In the following we also need the solution for the variance a qq , which is given by: 



<?<?' 

e~ 2Xt A 
a qq (t) = — {A„[(/i 2 + i/ 2 ) cosh 2vt + 2/iv sinh 2ut + to 2 } + -^|(cosh2z/t - 1) 
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H — A P9 (/x cosh 2ut + i/sinh2z/t — //)} + cr qq (oo), (17) 
where A qq = a qq (0) - a qq (oo), A pp = cr pp (0) - a pp (oo), A pq = cr pq (0) - a pq (oo) and 
M°o) = 2U , 2 1 2 2t(^ 2 (2A(A + /i) - uj 2 )D qq + D pp + 2m(A + p)D pq ), (18) 

ZTfl AyA — Ul — fl ) 

= 1 -{{mujfoj 2 D qq + (2A(A - n) - to 2 )D pp + 2mcu 2 (\ - (j,)D pq ),(19) 

ZAyA — UJ — fJ, ) 

a pq (oo) = ^——-^—— ((X + fi)(muj) 2 D qq + (A - fi)D pp + 2m{\ 2 - fi 2 )D pq ).{20) 

Note that in the case A > v, a q (t — > oo) = cr p (t — > oo) = and a qq {t — ► oo) = a qq (oo). 
In the case A < u, a q {t — > oo), a p (t — > oo) — > ±oo and cr w (t — > oo) = oo. 

3 Tunneling through an inverted parabola 

In order to calculate the tunneling probability through the inverted harmonic potential 
dHJ), we assume that initially the wave function of the system is a Gaussian wave packet 
centered at the left of the peak of the potential at q = 0, cr q (0) < 0, with a momentum 
cjp(0) > towards the potential barrier peak: 

^ = w^mw exv[ -^M {q - ff ' (0))2 + r- m - (21) 

Then the corresponding initial probability density is given by: 

" fet = °> = (2^(0))V3 ^1-2^(0)^ " ff »(°)» 2 l' < 22 » 
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Like in [3Q]-[32J, we can transform the master equation <^ for the density operator 
of a particle moving in the potential ( fl4l) of an inverted parabola into the following 
Fokker-Planck equation satisfied by the Wigner distribution function W(q,p, t) : 

dW v dW o dW d d 

sr = -sw " mw '*w + (A - "W qW) + (A + "W pW) 

d 2 W (PW_ 2D d 2 W 

dq 2 pp dp 2 pq dpdq 

For an initial Gaussian Wigner function, the solution of Eq. (1231) is 

1 



+D «nz- + D »i&r + 2D p^- ( 23 ) 



W(q,p,t) 



2nJa{t) 



x exp{-^y[a pp (t)(g - a q (t)) 2 + a qq (t)(p - a p (t)) 2 - 2a pq (t)(q - a q (t))(p - <r p (t))]},(24) 

which represents the most general mixed squeezed states of Gaussian form. Here 
cr q (t), cr p (t) and cr qq (t), cr pp (t) , a pq (t) are the expectation values and, respectively, the 
variances corresponding to the inverted parabola as given partly in Eqs. (fT5l) - ffTTj) and 

a(t) = a qq (t)a pp (t) - cr pq (t) 2 . (25) 

Since the dynamics is quadratic, then according to known general results, the initial 
Wigner function remains Gaussian. The density matrix can be obtained by the inverse 
Fourier transform of the Wigner function: 

f i q ~t~ q r 

< q\pW >= J dpexp(-p(q - q'))W(^—,p,t). (26) 

Using Eq. ( 124"1) . we get for the density matrix the following time evolution: 

2-KO qq (t) 2 &qq{t) 2 

~lL(a pp (t) - °^){q - q') 2 + p^-J^- ~ a q (t))(q - q') + U(t)(q - q')].(27) 
2k Vgqit) na qq (t) 2 h 

The initial Gaussian density matrix also remains Gaussian, centered around the clas- 
sical path, i. e. cr q (t) and cr p (t) give the average time-dependent location of the system 
along its trajectory in phase space. The wave function starts as a Glauber wave packet 
at t = on the left-hand side of the barrier and evolves as a mixed squeezed state at a 
later time. By putting q' = q in Eq. fT2"Tj) . we obtain the following probability density 
of finding the particle in the position q at the moment t: 

^ = {2,alm' 2 eM '^) {q ' aq{m (28) 



6 



This is a Gaussian distribution centered at cr q {t), which describes the classical trajec- 
tory of a particle initially at cr q (0), with initial momentum cr p (0) and variance cr qq (t). 

Using Eq. ( l28l) . the probability for the particle to pass to the right of position q 
at time t is given by 

P(Q, t) = p(q', t)dq' = f° 1 exp(-^-^^)rfg'. (29) 

We define the tunneling probability P(t) as the probability for the particle to be at 
the right of the peak at q = (beyond the barrier top): P(t) = P{q = 0,t). We obtain 

P(t) = -4= r ™ e~ u2 du = -(1 - erf( )), (30) 

where erf(x) is the error function with erf(x) = — erf(— x) and erf(oo) = l. 

From Eq. fl30|) we see that the probability P(t) depends only upon the classical 
motion of the average value of coordinate (wave packet center) and the spreading of 
the wave packet in the direction of the barrier. The final tunneling probability (barrier 
penetrability) is given by taking the limit t — > oo in P(t). In the present calculations 
we ignore the fact that a part of the wave packet has already tunneled through the 
barrier at t — 0. In general, this probability has a negligible value, but, in principle, 
in order to find the net penetration probability, it should be subtracted from the 
tunneling probability at time t. 

4 Evaluation of the penetration probability and 
analysis in dimensionless variables 



We will show that since u q {t) and Ja qq (t) are both proportional to the same exponen- 
tial factor as time approaches infinity, their ratio in Eq. (!30|) approaches a finite limit, 
which determines the final tunneling probability. Indeed, as t — > oo, we see from Eqs. 
( TrS]) and (|T7|) that cr q (t) and cr qq {t) behave like 

e -(A-i/)t 

^ < 31) 



and 



p -(A-i/)t 



° qq (t) ^ (~2^-) A + ^(°°)' ( 32 ) 
where we have denoted 

5 = m(/i + v)<t 9 (0) + a p {0) (33) 
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and 



A = m 2 (fi + v) 2 A qq + A pp + 2m(/i + u)A pq . (34) 
Then we obtain the following finite limit as t — > oo : 

- 6 (35) 



0~aa VA 
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if A < v (and the limit if A > v) and, therefore, the expression (I30I) leads to the final 
penetration probability [P = P(t — ► oo)) 

P = ^(l-erf(--^)) (36) 

if A < z/ and P = 1/2 if A > v. In the case A > u, u q (t) — > if t — * oo, that is, 
the system is located around the barrier, a qq {t) tends to a finite value a qq (oo) for any 
initial kinetic energy and in this case P(t) = 1/2. Let us consider the other case A < v. 
For S = the trajectory tends to the top of the potential barrier and if S is different 
from 0, then <J q {t) tends to oo or — oo. The trajectory which starts on the left-hand 
side of the barrier (cr g (0) < 0) with a positive initial momentum cr p (0) > will stay 
on the same side for 5 < (and then a q (oo) — > — oo) and will cross the barrier for 
5 > 0, i. e. if the initial kinetic energy allows to overcome the barrier (and then 
cT g (oo) — > oo). For a general fi and any A < u, the particle crosses the barrier when 
<7 P (0) > —m(fi + u)cr q (0) and will stay on the same side when c p (0) < — m(fi + v)a q (0). 
The barrier penetrability is larger than 1/2 if the particle classically can overcome the 
top of the barrier, it is smaller than 1/2 if the particle cannot cross the barrier and it 
tends to 1/2 if the position uncertainty a qq is very large. 
For /i = 0, the values of A and 5 become {v = uj): 

5 = mwa q {0) + a p {0) (37) 

and 

A = m 2 u 2 (T qq (0) + appiO) + 2mwa m (0) qq pp ^. (38) 

A — uj 

The particle crosses the barrier if cr p (0) > -mwfr 9 (0) and A < uj. In this case, if 
A increases, then the ratio 5o/(2A ) 1 ^ 2 and the penetration probability P decreases. 
This means that if dissipation increases, then the probability P decreases up to a value 
of 1/2. If (Jp(0) < — muj<y q (0), the particle can not cross the barrier. In this case the 
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penetration probability P increases with the dissipation A up to 1/2. At A > uo the 
wave packet sticks in the barrier region. 

We now introduce the dimensionless variables: z the scaled initial position, v the 
scaled initial momentum, e the scaled dissipation coefficient and r the scaled inverse 
wave packet size, defined as follows: 



z= " g(0) , v= 3% c=A r = \ 2m ". (39) 
/MO) ma;a,(0) w ^ ( ) 



With these notations and considering a thermal bath modeled by the coefficients of 
the form (jSJ), the penetration probability takes the following form for the case fi = 0: 

- ^ - A» ■ 5 (1 - erf( - T/ rft „^ )) - (40) 

V2 J 1 + r 4 r 2 coth — — 

V e - 1 2£;T 

We took into account that 

M0)M0) = t?/4, (o- M (0) = 0). (41) 

If n = 0, then we have < e < 1. 

For /i ^ 0, in the case of a thermal bath, the expression (j34p takes the form 
(a p? (0) = 0) 

A = m 2 (/i + i/)V w (0) + a pp (0) 

km . , . u; 2 [A 2 + A(u + 1/) + uz/1 . , ?ia; , ,„ s 

{M/* + *0 + 1 x9 , V-^} coth——. (42) 



oj Lf ^ ' \ 2 -iu 2 -fi 2 J 2kT 

With the notations (j3"9"l) and introducing also the notation 



7 = ^, (43) 



the penetration probability takes the following form: 



^o = ^(l-erf(--J=)), (44) 



with 

5 2(7 + v / l + 7 2 + 



A [(7 + v^) 2 + H - 2( 7 + VTTf) (£ 7 M±f + V coth ^ 



.(45) 
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If \i 7^ 0, the inequalities A > /x (see Eq. (JSJ) and A < v lead to the following restrictions 
on the dimensionless variables: 



7 < e < y'I + 7 2 . (46) 
The initial energy of the particle associated with the Gaussian wave packet (21) is 

E=<H> \ t=0 = ^ pp (0) —a qq (0) + —a 2 p (0)-—a 2 q (0) +^(0)^(0) (47) 

and in terms of the dimensionless variables (39) it looks 

£=^[r 4 -l + *V-l)] + Y7^ ( 4 8) 

If ^ < it is a sub-barrier initial energy and if E > it is an energy above the 
barrier. In terms of the same dimensionless variables, the condition that a classical 
particle does not have enough initial kinetic energy to pass the potential barrier can 
be written: 

v > -1, if p. = (49) 

and 



v> -(7+ y/l + 7 2 ), if// t^O. (50) 

With these two conditions and by taking < r < 1 (which assures that the initial 
fluctuation energy is negative), the total initial energy E (47) is always negative. This 
corresponds to the case of the sub-barrier energy, relevant to the quantum tunneling 
problem. The examples provided in the following figures reflect just this situation. 

For /i = 0, Figs. 1 and 2 show the dependence of the tunneling probability on the 
scaled dissipation e and the temperature T of the thermal bath, via coth(huj / 2kT) , for 
fixed values of the scaled initial position z, momentum v and wave packet size r. 

In the next four figures, we consider /i ^ 0. Figs. 3 and 4 show the dependence 
of the penetration probability on the scaled dissipation and on the parameter 7 for a 
fixed scaled initial position z, momentum v and wave packet size r at the temperature 
T = 0. 

In Figs. 5 and 6 we give the dependence of the penetration probability on the 
scaled dissipation and temperature, at fixed values of z, v, r and 7. The presented 
dependence of the penetration probability on these variables can be summarized in 
the following conclusions: 

1) When the scaled initial momentum \v\ is increasing, then P is increasing up to 
1/2 if the particle does not have enough kinetic energy to pass the potential barrier. 
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Figure 1: Dependence of tunneling probability P on the scaled dissipation e = \/u 
and the temperature T of the thermal bath, via coth(huj/2kT), for /i = and for fixed 
values of the scaled initial position z = —3, scaled initial momentum v = —0.5 and 
scaled inverse wave packet size r = 0.5. 
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Figure 3: Dependence of tunneling probability P on the scaled dissipation e and pa- 
rameter 7 = fi/u, for the temperature T = of the thermal bath and for fixed values 
of the scaled initial position z = —3, scaled initial momentum v = —0.5 and scaled 
inverse wave packet size r = 0.3. 
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epsilon 8.015 

8 . 02 u 

Figure 5: Dependence of tunneling probability P on the scaled dissipation e and the 
temperature T of the thermal bath, via coth(hu/2kT), for 7 = 7.99 and for fixed 
values of the scaled initial position z = —3, scaled initial momentum v = —0.5 and 
scaled inverse wave packet size r = 0.5. 




1 . 1 

Figure 6: Same as in Fig. 5, but with 7 = 0.97, z = —9, v = —0.9 and r = 0.5. 
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The same conclusion is valid for the variable r, i. e. if the initial width of the Gaussian 
packet is decreasing, then the penetration probability is increasing. 

2) If the scaled initial position \z\ is increasing, then P is decreasing from 1/2 to 
if the particle does not have enough kinetic energy to pass the potential barrier. 

3) The penetration probability is increasing from to 1/2 with dissipation and 
with coth.(hu/2kT) and, therefore, with the temperature. For the case /i ^ 0, the 
probability P is decreasing with \x. 

In conclusion, the dependence of the tunneling probability on dissipation is not 
simple. When the particle does not have enough kinetic energy to pass the parabolic 
barrier, which is the relevant case to the quantum tunneling problem, the dissipation 
enhances tunneling. 

5 Summary 

In the framework of the Lindblad theory for open quantum systems, we have formu- 
lated the motion and the spreading of Gaussian wave packets in an inverted oscillator 
potential. We have obtained analytic solutions of evolution in time of the wave packets 
and of the barrier penetrability. Since the wave packets spread in time according to 
the same law of evolution as their center moves, the value of barrier penetrability is in 
general different from 1/2. The inverted oscillator potential has an important physical 
relevance, since it can constitute a guide how to treat more physically realistic poten- 
tials, like third order and double-well potentials (23, [33] or joined inverted parabola 
and harmonic oscillator potentials [34] , in order to be applied in nuclear fission and in 
molecular or solid state physics. 
Acknowledgments 

One of us (A. I.) is pleased to express his sincere gratitude for the hospitality 
at the Institut fiir Theoretische Physik in Giessen. A. I. also gratefully acknowledges 
financial support by the DA AD (Germany). 

References 

[1] A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 211 (1981) 
[2] A. O. Caldeira and A. J. Leggett, Ann. Phys. (N. Y.) 149, 374 (1983) 
[3] A. J. Bray and M. A. Moore, Phys. Rev. Lett. 49, 1545 (1982) 



14 



[4] S. Chakravarty, Phys. Rev. Lett. 49, 681 (1982) 

[5] G. Schdn and A. D. Zaikin, Phys. Rep. 198, 237 (1990) 

[6] S. Chakravarty and A. J. Leggett, Phys. Rev. Lett. 52, 5 (1984) 

[7] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg and W. 
Zweger, Rev. Mod. Phys. 59, 1 (1987) 

[8] U. Weiss, Quantum Dissipative Systems (World Scientific, Singapore, 1992) 

[9] H. Dekker, Phys. Rev. A 38, 6351 (1988) 

[10] G. W. Ford, J. T. Lewis and R. F. O' Connell, Phys. Lett. A 158, 367 (1991) 

[11] M. Razavi and A. Pimpale, Phys. Rep. 168, 305 (1988) 

[12] H. Hofmann, Phys. Rep. 284, 137 (1997) 

[13] A. Schmid, Ann. Phys. (N. Y.) 170, 333 (1986) 

[14] A. Widom and T. D. Clark, Phys. Rev. Lett. 48, 63 (1982) 

[15] R. Bruinsma and P. Bak, Phys. Rev. Lett. 56, 420 (1986) 

[16] A. J. Leggett, Phys. Rev. B 30, 1208 (1984) 

[17] M. Razavy, Phys. Rev. A 41, 6668 (1990) 

[18] K. Fujikawa, S. Iso, M. Sasaki and H. Suzuki, Phys. Rev. Lett. 68, 1093 (1992) 

[19] E. G. Harris, Phys. Rev. A 48, 995 (1993) 

[20] L. H. Yu, Phys. Rev. A 54, 3779 (1996) 

[21] J. Ankerhold, H. Grabert and G. L. Ingold, Phys. Rev. E 51, 4267 (1995) 

[22] J. Ankerhold and H. Grabert, Phys. Rev. E 52, 4704 (1995); 55, 1355 (1997) 

[23] E. B. Davies, Quantum Theory of Open Systems (Academic Press, New York, 
1976) 

[24] G. Lindblad, Commun. Math. Phys. 48, 119 (1976) 

[25] H. Spohn, Rev. Mod. Phys. 52, 569 (1980) 

15 



[26] E. Stefanescu, A. Sandulescu and W. Greiner, Int. J. Mod. Phys. E 2, 233 (1993) 

[27] G. G. Adamian, N. V. Antonenko and W. Scheid, Phys. Lett. A 244, 482 (1998) 

[28] A. Isar, A. Sandulescu, H. Scutaru, E. Stefanescu and W. Scheid, Int. J. Mod. 
Phys. E 3, 635 (1994) 

[29] A. Sandulescu and H. Scutaru, Ann. Phys. (N. Y.) 173, 277 (1987) 

[30] A. Isar, Helv. Phys. Acta 67, 436 (1994) 

[31] A. Isar, W. Scheid and A. Sandulescu, J. Math. Phys. 32, 2128 (1991) 

[32] A. Isar, A. Sandulescu and W. Scheid, Int. J. Mod. Phys. B 10, 2767 (1996) 

[33] G. G. Adamian, N. V. Antonenko and W. Scheid, Nucl. Phys. A 645, 376 (1999) 



[34] S. Misicu, quant-ph/9906012 



16 



